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Ordinary tensor models of rank D > 3 are dominated at large N by tree-like graphs, known as 
melonic triangulations. We here show that non-melonic contributions can be enhanced consistently, 
leading to different types of large N limits. We first study the most generic quartic model at D = 4, 
with maximally enhanced non-melonic interactions. The existence of the 1/N expansion is proved 
and we further characterize the dominant triangulations. This combinatorial analysis is then used to 
define a non-quartic, non-melonic class of models for which the large N free energy and the relevant 
expectations can be calculated explicitly. They are matched with random matrix models which 
contain multi-trace invariants in their potentials: they possess a branched polymer phase and a 2D 
quantum gravity phase, and a transition between them whose entropy exponent is positive. Finally, 
a non-perturbative analysis of the generic quartic model is performed, which proves analyticity in 
the coupling constants in cardioid domains. 
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I. INTRODUCTION 

An interesting way to define and work with continuous random geometries (for instance, the continuous random 
tree) is to start with random discrete spaces (for instance, trees) whose sizes can become arbitrarily large. In a suitable 
scaling limit, such random discrete objects can converge in a precise (Gromov-Hausdorff-Prokhorov [IJ[2j, or GHP for 
short) sense to a continuous random space of a certain universality class. 

It is standard for physicists to think of discrete random spaces as combinatorial objects generated by the Feynman 
expansion of some integrals over tensors. A tensor of rank D is a D-dimensional array of (complex) numbers. A 
position in such an array is labeled by a D-uple of integers each ranging from 1 to N. Of particular interest to us 
are random tensors equipped with a (/(Vj^-invariant distribution. If the Feynman expansion can be organized as a 
series in 1/N, then the large N limit exists and selects a class of combinatorial objects (the random spaces) to be 
summed over. 

A tensor of rank 1 is a vector and the large N limit of vector models consists of plane trees (in the intermediate 
held representation) whose GHP limit is the continuous random tree (CRT, for short) [3]. The CRT, also known as 
the branched polymer universality class, has Hausdorff dimension 2 and spectral dimension 4/3. 

The situation changes drastically for tensors of rank two, i.e. random matrices. Indeed, matrix models provide a 
quantization of two-dimensional gravity coupled to conformal matter gUj. Such matrix integrals typically generate 
maps (possibly decorated with matter) weighted by N 2 ~ 2g where g is the genus of the map [7j. This topological 
expansion is dominated at large N by planar maps (of spherical topology). As random spaces, planar maps (without 
critical matter) converge in the GHP sense to the Brownian sphere [H [9], a random space with Hausdorff dimension 
4 and (expected) spectral dimension 2. 

In higher dimensions matrix models generalize to tensor models mm and group held theories H2HH]- However 
such models did not admit a large N limit. The situation changed with the advent of colored US; HJ and U(N) D - 
invariant m tensor models. They admit a 1/N expansion [151 - 121] . This perturbative Feynman expansion furthermore 
sums over all manifolds of dimension D and a restricted class of pseudo-manifolds [22]. The weights of those objects are 
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moreover simple discretizations of the Einstein-Hilbert action |23l . coinciding with Euclidean dynamical triangulations. 
Hence tensor models seem relevant for the quantization of gravity in dimensions higher than 2 |24fl26j . They can 
support matter fields [2ZH33 and have statistical mechanics applications m- 

The first analytic solution of tensor models at large N, for D > 3, matched the expectations from the numerical 
analysis of Euclidean dynamical triangulations, [131 [32] . Mainly, they are dominated by melonic graphs which are 
representatives of the universality class of branched polymers [53] . The coupling to matter reproduces the critical 
exponents of multi-critical branched polymers. Notice that the colored and invariant families of tensor models admit 
double-scaling limits for D < 6 [34fj37l . which are also made of trees. We would like to point out that trees can serve 
as building blocks for more elaborated objects. A key idea is to decorate them cleverly with a certain density of loops, 
for instance by additional data such as Schaeffer’s “well-labelings” [35] (the process which leads to the Brownian 
sphere), or by mating two trees (a process which leads to a kind of Brownian sphere decorated with an SLE called 
a “peanosphere” [39]). The present article precisely shows a way for tensor models to go from trees to planar maps 
through a phase transition which exhibits a proliferation of planar maps along a tree. 

There is indeed more to tensor models than melonic graphs. For instance, one can trivially build a square matrix 
from a tensor of rank 4 and equip it with an ordinary matrix model, which would converge at large N to the Brownian 
sphere. Importantly, this requires the distribution on the random variables to have a different behavior with N than 
those dominated by melonic graphs. Furthermore, the boundary observables of tensor models are far more richer 
than those of matrix models: instead of collections of loops, they are labeled by /3-regular graphs. Most of them 
usually contribute to sub-leading corrections, even when they explicitly enter the measure. However, it might be 
possible to modify the behavior of the measure with N so as to enhance such non-melonic observables. This idea was 
first developed in [40 ]. It showecQthat while some models with non-melonic distributions simply extend the range of 
applicability of the universality theorem of [41] , some other models can escape it and converge to some random space 
that is not the CRT. 

It was in fact noticed in the last section of |41| that every tensor observable which appears in the measure should 
admit a range of admissible scalings with N for which the 1 /N expansion exists. For a melonic invariant, this range is 
restricted to a single possible value, which is the one used in [mill]. But non-melonic invariants can lead to different 
types of 1/N expansion, as established and exploited in ■ The present article pushes this approach further: we start 
the analysis of a new kind of 1/iV expansion for tensor models. It is based on melonic interactions, and non-melonic 
ones which in contrast to the existing literature (except for [3D]) are appropriately rescaled to their largest admissible 
scaling. We say that those observables are enhanced and call such models maximally rescaled. 

In Section [XT] we explain the notion of maximally rescaled observable and briefly review basic facts about tensor 
models (including the structure of the Feynman graphs and a comparison of the present approach with 140 ] ). 

In SectionJTTT] the most generic quartic, maximally rescaled model at D = 4 is analyzed. The combinatorial analysis 
is performed thanks to the intermediate field method. It proves the existence of the 1/iV expansion and characterizes 
the graphs which dominate the free energy at large N. 

In Section [Tvl we introduce another generic class of maximally rescaled models whose measures involve non-melonic, 
non-quartic invariants which we call trees of necklaces (still at D = 4). The motivation is that these observables 
are boundary observables of the dominant graphs of one of the quartic models. This makes possible to use the 
combinatorial analysis of Section |III| to identify the graphs which dominate at large N. Combined with Schwinger- 
Dyson equations (or Tutte equations), that leads to their enumeration which allows the large N evaluation of the 
free energy (in particular its critical entropy exponents) and expectations of observables. While the critical behaviors 
which arise in this class of models are known from the literature on matrix models with multi-trace invariants [42H46] , 
our analysis paves the way for the future analysis of more complicated maximally rescaled models which would better 
take advantage of the richness of tensor observables. 

Finally, Section |V] offers a non-perturbative analysis of the quartic maximally-rescaled models for any rank D. It 
establishes uniform analyticity and Borel summability at large N in cardioid domains of the couplings. That is a 
non-trivial extension of the results of [47] which establishes similar results for quartic models of any rank, but without 
enhancement of non-melonic observables. 


This analysis was also required in order to define tensor models for “rectangular” tensors. 


l 
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II. MODELS WITH MAXIMALLY ENHANCED INTERACTIONS 
A. Maximally enhanced observables 

Consider a Hermitian inner product space V of dimension N and {e n \n = 1,... ,N} an orthonormal basis in V. 
The inner product allows us to identify the dual of V and its complex conjugate V. hence not to distinguish between 
{e n \n = 1,..., N} and its dual basis, nor between upper and lower indices. A complex tensor of rank D is a multilinear 
form T : V® D —► C, which can be written in components as 

T = 'y ' T ni ...n D e ni 0 * * * 0 • (1) 

It is important that T ni ... no has no symmetry properties. Hence its indices have a well defined position. We call the 
position of an index its color. We have therefore a set T> of D colors. From now on we specialize to rank D = 4 hence 
tol> = {l,... 4}, and T = J2 ni ,n 4 Tm-m e ni 0 • • • 0 e„ 4 . The dual tensor T is defined by 

T= T m-ru e ni 0---0e n4 . (2) 

n 1 ,...,n 4 

The unique quadratic invariant is the (scalar) Hermitian pairing of T and T which writes: 

T ~d T = ^ T ni ... ni T ni ... n4 . (3) 

7li,n2,n3,n4 

Observables are rnonic, homogeneous polynomials invariant under U(N) D . Given a finite set of such observables, 
{Hj(T, T)} ieJ , we define the action S, the partition function Z and the free energy F as 

S( T,T) = T -a T + ^2uBi( T,T), 
iei 

- ( 4 ) 

exp(—F 1 ) = Z = f exp(-iV Q A(T,T)). 

J n 

The exponent a which parametrizes the measure is usually fixed by the following requirements: 

• the free energy is bounded, for any choice of {Bi} ieI , by a polynomial bound at large N of type |F| < N K , 

• there exist interactions so that expanding the free energy perturbatively in the coupling constants 

around the Gaussian measure, there are infinitely many Feynman graphs with the same exponent of N, 

• the coupling constants are independent of N. 

These criteria uniquely select a = D — 1. The first two requirements are crucial. The first one ensures that the free 
energy, appropriately rescaled like F/N K , is not divergent at large N. This requires a not to be too small. The 
second one is necessary to avoid trivial models where only a finite number of Feynman graphs would contribute at a 
fixed order of the 1/N expansion. This requires a not to be too large. 

However, the last requirement does not seem equally important. In fact, it can be relaxed to the benefit of getting 
new large N behaviors, as we will show in the present article. The situation is slightly similar to the case of matrix 
models. There D = 2 and a = 1 for actions which are sums of single-trace invariants. Multi-trace invariants can be 
included provided they come with TV-dependent weights, otherwise the free energy is divergent. 

To see how that works, we need to give a bit more information on the usual large N limit of tensor models. The free 
energy goes like N D and it is dominated by the melonic observables, i.e. F only depends at large N on the coupling 
constants of melonic polynomials. In other words, the second criterion above is only satisfied for melonic interactions. 

If B is non-melonic, the H-action Sb = T-jiT + ts B( T, T) leads to a free energy Fb which is independent of 
tB at large N. However, there may exist an IV-dependent weight ts = N aB ts , with ts A-independent, such that 
the associated H-action preserves the boundedness of the free energy while leading to a non-trivial (and non-melonic) 
model. If such an as >0 exists, we say that B can be maximally enhanced. Here “maximally” refers to the fact that 
if the value of as is increased, boundedness of the free energy is lost. 

It is quite simple to see that such models do exist. For instance, set D = 4 and consider T as a matrix of size 
N 2 x N 2 between pairs of tensor indices. Then, matrix model building applies as usual (and leads to a 2D quantum 
gravity phase at large N). 
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FIG. 1. The seven invariants at rank 3 and order 6 : there are two different melonic families each with three members (shown 
left and center) and the non-planar bipartite complete graph A 3,3 (shown right). 


This was already observed in [3D], with a different but equivalent strategy. Instead of keeping a = D — 1 fixed and 
enhancing non-melonic polynomials, it was chosen to decrease a < D — 1 and simultaneously add some TV-dependent 
weights on melonic observables so that F does not diverge. It is easy to check that upon a rescaling of T, this is 
equivalent to the approach we have just described. 

However, our present work improves [401 on many aspects. For instance, the models described in [3D] are labeled 
by “slices of colors”, i.e. a choice of partition of the D indices of the random tensor. As a consequence, observables 
which only differ by a permutation of some tensor indices would receive in |40l different TV-dependent weights. There 
is no such partition here and all observables which differ by a relabeling of the indices receive the same enhancing. 

Furthermore, in the context of the color slices, it was proved in [40] that the large TV limit is always of the same type 
as in tensor models with their standard scaling, except if the action contains some generalized matrix-like observabletrl 
However, the large TV limit of the latter models was not studied, and neither were models containing both melonic 
and matrix-like interactions. In the present article, we will give a much more comprehensive analysis of what happens 
when melonic observables are mixed with certain matrix-like observables. This mixing means that those observables 
will all appear in Feynman graphs. We will then see that this can be used to include new observables which result 
from gluing melonic and matrix-like observables together in the action of the model. 


B. Feynman Graphs 

The ordinary perturbative expansion is obtained by expanding the exponential of the non-quadratic terms of S in 
Q (called the interactions of the theory) into a series in the couplings then (illegally^ commuting the sum 

of the series with the remaining Gaussian integral. The resulting series is indexed by Wick contractions, also called 
Feynman graphs which we now describe. 

First, we recall that invariant polynomials can be represented as bipartite, D-regular graphs with edges colored 
in {1,... ,D} such that no two edges incident to the same vertex carry the same color. A white (or black) vertex 
represents T (or T) and an edge of color c between two of them represents the contraction of their indices in position 
c. We call such graphs bubbles. Some bubbles are on display in Figure [l] 

The Feynman graphs generated by tensor integrals consist of bubbles chosen in arbitrary (finite) numbers among 
those which represent the invariant polynomials connected via additional edges to which the color 0 is 

attributed. The resulting graphs are bipartite, ( D + l)-regular graphs whose edges have a color in {0,..., D} such 
that the edges incident to a given vertex all have different colors. We call them vacuum Feynman graphs and the free 
energy expands as a sum over the connected ones. 

If Q is a connected vacuum Feynman graph generated by the logarithm of the integral 0. containing rii(Q) bubbles 
of type Bi, for each i G /, then the amplitude associated to Q reads 

A(Q) = N F( - s) - aE °(S)+«£ ieJ ni(S) J|(-i i ) n< ( £! ), (5) 

iei 

which also works in the maximally rescaled case where the couplings t t may have an TV-dependence. Here Eq(Q) is 
the number of edges of color 0 (it equals the number of white vertices of Q). The quantity F(Q) is the total number 
of faces. A face of colors (0c), for c £ {1,..., D}, is a cycle whose edges have the colors 0 and c only. One has the 
partition F{Q) = J2c=i F Oc(G)- 


2 These observables are defined technically as having a single face of colors (i,j) and a single face of colors ( k,l ), where {i,j, k,l} = 
, { 1 , 2 , 3 , 4 }. 

3 The legal version of this procedure is constructive theory, see Section Ivl 
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FIG. 2. The quartic invariants at rank 4: there are four types of melonic graphs (one of them shown left, the others being 
obtained by permuting the colors) and three types of necklaces (one of them shown right, the others being obtained by color 
permutations). 

III. ENHANCED QUARTIC MODELS 

At rank 3, the three quartic (two black and two white vertices), connected bubbles are all melonic (and the associated 
model, restricted to a single coupling, i.e. invariant under permutations of the colors, was built non-perturbatively in 
(48)). The first non-trivial maximally rescaled tensor model with quartic interactions has to be looked for at D = 4. 
At that rank there are seven quartic bubbles, four of which are melonic and three non-melonic, as depicted in Figure 
[21 We call the non-melonic bubbles necklaces. 


A. Rank 4 Quartic Tensor Models 


We simplify slightly the notations of [47] , to adapt to the rank 4 case. The normalized Gaussian measure associated 
to the unique quadratic invariant © is 


duo= (n 


dT n dT n \ _jv a T- D T 
2m 


( 6 ) 


A connected quartic t/(iV) 4 -invariant polynomial Be is specified by a partition of V = {1,2,3,4} into two non¬ 
trivial subsets C CD and its complement D\C. The associated invariant can be written as a trace over the indices 
of C of a partial scalar product over the indices of D \ C or vice-versa [47] 


B C (T,T) = Tr c 


[T -D\c Tj 


•c [T " d\c T| 


= Tr-p_c 


[T - c T] . V \ C [T - c T] 


(7) 


where we denoted x the product of operators from U® c to V® c . 

There are therefore 7 quartic connected invariants at rank 4, the 4 melonic ones which correspond to type 1-3 
partitions in which one subset is a singleton (Figure [2] left hand side), so to Bc t with C t = ji}, i = 1,2, 3,4 and the 
three non-melonic, necklace partitions (Figure 2] right hand side) which correspond to type 2-2 partitions, hence to 
Bc u with Cu = {1,*} with i = 2,3,4. Remark that such necklace interactions can be interpreted also as ordinary 
matrix traces on an N 2 space, namely V <§) V. 

More precisely in components we have: 


BcA T,T) = 


E 


T T / / / T / / / / T / 

nin 2 n 3 n 4 - 1 - 71^712^3714 




and three similar formulae for Bc 2 , Bq 3 and Bc 4 , replacing the special position 1 by 2, 3 and 4. Also 


Bc 12 ( T,T)= £ 


D mn2n 3 n 4 D nin2n '^ n '^T n ^ n ' 2n '^ n ^T n ^ n '^ n3n4 


( 8 ) 


(9) 


Tl i,... ,724,71^ ,...,714 


and two similar formulae for Bc 13 and Bc 14 replacing the special pair 12 by 13 and 14. 

The standard general quartic tensor model at rank 4 is then the (invariant) perturbed Gaussian measure 


d» standard = d„ 0 e- N3 ^KB Ci CT) -N 3 EUb iBcij (T,T) . 


( 10 ) 


The standard melonic color-symmetric quartic model at rank 4 corresponds to Ai = A 2 = A 3 = A 4 = g\ and to 
A 12 = A 13 = A 12 = 0, hence it has only one coupling constant g\. It has been proved Borel summable uniformly 
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in N for g\ in a cardioid domain in [45]. The standard color-symmetric quartic model at rank 4 corresponds to 
Ai = A 2 = A 3 = A 4 = g 1 , and to A 12 = A 13 = A 12 = 52 , hence it has two coupling constants. It has been proved Borel 
summable uniformly in N for gi and g% in cardioid domains in m- 


The maximally rescaled general quartic model at rank 4 corresponds to changing in (10) the scaling of the necklace 
terms from N 3 to N 4 . This corresponds to their more natural scaling as matrix traces over some V (g> V spacc^J and 
is also the maximal rescaling for which the 1/7V expansion of the tensor model still exists, 


d^max — dl>i) & 


-N 3 Et= 1 EB Ci (T,T) -N 4 Et = 2 Ai i B c u (T.T) 


( 11 ) 


The maximally rescaled restricted quartic model at rank 4 (in short called from now on the restricted model ) is the 
same model but without the interactions Bc 13 and Bc Xi ' 


dg res = dg 0 e -^Eti^(T,T) -at 4 a 12 b Ci 2 (t,t) _ 


( 12 ) 


We are interested in the generating function of the moments and the cumulants k of the measure d/i, which are 
defined as: 


Z(jJ) 



E n T n J n +E n TnJn 


n(T ni T ni ...T np T np ) 


(in Z(J,J)^j 

(irJj, ! ()rJ n . . .<)J rip ()J rlp 


J=J—0 


(13) 


In order to compute the cumulants of g we need to compute the logarithm of the generating function Z(J, J). This 
is where the field-theoretic notion of Feynman graphs is useful: the function Z(J,J) expands as a power series in 
g 1 and 32 indexed by Feynman graphs, and its logarithm expands as power series in gi and <?2 indexed by the same 
Feynman graphs but connected. 


B. Intermediate-field representation 



FIG. 3. Intermediate Field Graphs 


4 This can be seen by changing T with T /y/N. Both the Gaussian term and the necklace interactions then get a global factor IV 2 , 
indicating a matrix model over V <S> V where V has dimension TV 2 . More details are provided in the following section. 
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Quartic models are most conveniently studied in the intermediate field (IF) representation (see Figure [3]). In our 
case it has been defined in detail in HU EH], so that we just give an informal description sufficient for our present 
purpose. The graphs of this IF-representation are in fact maps , i.e. there is a cyclic order of the edges incident to 
every vertex. Each IF vertex corresponds to a former cycle of 0-edges in the ordinary or O-representation, hence 
should be represented as a 4-stranded loop, with strands of colors 1,2 3 4 corresponding to the position of the tensor 
indices. The corners (between two incident edges) of the IF representation represent the former edges of color 0 
of the O-representation. Bubbles of the type Be- of the O-representation become monocolored edges labeled with 
the color f, opening the single strand of color f of the IF vertex. Bubbles of the type Bc u in the O-representation 
become bicolored edges of color type 1 f, opening simultaneously both strands of colors 1 and i = 2,3,4 of the IF 
vertices. Hence an IF graph of the full quartic model is a map with edges labeled by seven different color types (four 
monocolored possibilities and three bicolored types 1 i) and an IF graph of the restricted quartic model is a map with 
edges of five different labels (four monocolored or the pair 12). 

Finally external pairs of half-edges become cilia on the vertices in the IF representation; the exact correspondence 
is quite subtle and described in detail in [47] I48j, to which we refer the reader; let us simply remark that a cyclic IF 
vertex can bear either 0 or 1 ciliurn but no more. 

Also rather than computinga general cumulant, we may be interested in the expectation value of a particular (gen¬ 
erally not connected) invaranlr] The relationship between cumulants and invariants is subtle and involves Weingarten 
functions (47] EH] • For simplicity we shall consider most of the time the simples^ function of the theory, namely the 
2 -point function 


G 2 


(in Z(J, J)) 
BJ 8 7 


(14) 


It corresponds to the expectation value of the unique invariant of order p = 1 defined in (|3j) . It expands over IF 
connected ciliated maps with a single ciliated vertex which provides a natural root for the entire map, and Weingarten 
factors are trivial. More complicated observables can be also treated by the same method. 

Consider a map A4 of the IF representation associated to a graph Q of the O-representation, and let us translate 
formula ([5| in the IF language for the full quartic model We denote Ei, i = 1,2,3,4, the number of edges in 
A4 with color f, and Eu the number of bicolored edges of color type (If), i = 2,3,4. Moreover, the faces of color Of 
in Q become the faces of the submaps Ai; which contain all the edges of color i and color type If and only them (a 
vertex with no such incident edge is an isolated vertex in A4; which has one face). We can still think of them as faces 
of colors Of since the color 0 labels the corners of A4. Denote F a , the number of faces of colors Of, then 


4 4 

a(m)= n(-^ N~ n(M ^ 

i =1 1=2 

where the exponent fl(AI) of 1/A in the amplitude for a map A4 is 

4 4 4 

Q(M) =3j2Ei + 2j^E li ~J2 F oi- 

2—1 C—Z 2=1 


(15) 


(16) 


C. The 1/N expansion 

Before moving on to the main results of this section, we prove the following lemma which will be used throughout the 
section. We recall that a cut-edge , or a bridge, is an edge of a connected map whose deletion transforms a connected 
map into two connected maps. 

Lemma 1. Let A4 be a connected map and e £ A4 an edge that is not a cut-edge. Denote A4 e the connected map 
obtained by removing e. Then, 

1. either e is a monocolored edge and H(A4) > fl(Af e ) + 2, 


5 Tensor invariants are both the interactions and the observables of a tensor model. 

6 The function log which expands into connected maps without cilia (p = 0), is apparently even simpler, but in fact the lack of 

a natural root slightly complicates its combinatorics. Also it requires to factorize first an overall factor A 4 in each vacuum map before 
attempting to compute the 1/A expansion. This A 4 factor is the analog of the usual volume to factor in a thermodynamic limit and 
corresponds to the computation of an intensive rather than extensive quantity. 







2. or e is a bicolored edge and Q(M) > Q(M e )- 

Moreover, if instead e € Ml is a monocolored cut-edge which separates Mi and M 2 , then fl(M) = fl(.Mi)-|-f 2 (A 42 )+ 4 . 

Proof Denote SO = 0(M e ) — O(M) for any quantity O defined on M e and M. We first prove that if the color 
label of e contains the color i then |<5Foi| < 1- Indeed, there are at most two faces of color 0* passing along e in M. If 
there are two distinct faces, they are merged together upon deleting e. If there is a single face along e, the deletion 
can at most split it into two. 

• Then, if e is monocolored of color i, Sfl = 3SEi — SF^i. Therefore Sfl < 38Ei + |<5Toi|- With SEi = —1, this 
implies Sfl < — 2 . 

• If e is bicolored with colors 1 i, then SLl = 2 SEu — SF 0 i — 5F oi < 2 SEu + |<5W 0 i | + \SF oi \ = — 2 + 1 + 1 = 0. 

If e is monocolored of color i and separates Mi and M 2 , then removing e increases F oi by 1 and decreases E, by 
1, all other quantities being unchanged. Therefore SLl = —SF oi + 3SEi = —4 and Ll(M e ) = S2(A4i) + fl(M 2 ). □ 

Establishing the existence of the 1/IV expansion amounts to proving that Q(M) is positive for all maps, up to an 
overall IV 4 for vacuum maps. 

Theorem 1 . For any vacuum connected map Ll(M) > —4; for any two-point connected map Q(M) > 0 and for any 
connected map for a higher order cumulant fl(M) > 2 . Moreover, for any vacuum tree fl(T) = —4; for any 2-point 
tree fi(T) = 0 . 

The same theorem has been already established (e.g. in [47] ) for the standard scaling fl'(M) = 3(Y)f =1 Ei + 
E\-,) — Eq,;. But this standard scaling is obviously strictly greater than Ll(M) if Y^c=iEu > 0. Hence 

Theorem [l] is a non-trivial new result. The fact that any cumulant of order higher than 2 is suppressed as IV —> 00 is 
similar to the standard asymptotic freeness in matrix or tensor models [41] . 

Proof We first show that for any closed connected map M, and any spanning tree T C A4, Ll(M) > S2(T). Then, 
we show that Q is constant on trees. 

Let T be a spanning tree in M . We delete the edges in the complement and track the variations of H in the process 
using Lemma [lj Order the edges in M \ T arbitrarily from 1 to S = \M \ T| and define the finite sequence of maps 
(A4( s )) s= o,...,s with M ^ = M and for s = 1 ,..., S, M^ as M ^ s_1 ^ with the edge e s removed (and M ^ = T). 
None of the edges in M \ T are cut-edges, so that all the maps M^ s) are connected. Moreover, Lemma [l] applies at 
each step. If e s is monocolored fl(M^) + 2 < H(A4^ S_1 ^), and if e s is bicolored, f2(A4^- ) ) < H(A4^ -1 ^). Clearly this 
leads to H(T) < H(A4). 

Notice that the bound obtained upon a monocolored deletion is strict, fl(A4C- ) ) < f2(A4^ s_1 ' ) ), while it is soft for 
bicolored deletions. This reflects the fact that in the strictly monocolored case, i.e. the melonic case, only trees can 
dominate the large N limit, while the bicolored edges may (and actually do) change that result. 

Now we prove that all vacuum trees scale the same, then get a bound for trees with a fixed number of cilia. 

• If M is the vacuum map reduced to a single vertex, it has amplitude ./V 4 (one face for each submap Mi of color 
i). Adding a monocolored edge of type i will add one vertex, hence three new faces (those for j 7 ^ i ) and a factor 
N ~ 3 , hence is neutral in N. The same thing holds for adding a bicolored edge of type 1 i: it adds two new faces 
and a factor N~ 2 , hence is neutral in N. Then, by induction, every vacuum map that is a tree scales like IV 4 . 

• Every 2-point tree map is obtained by adding a single ciliurn to a vacuum tree; this ciliurn opens exactly 4 
strands, hence faces, hence every 2-point map that is a tree scales at 0(1). Finally a tree for any cumulant of 
higher order is obtained by adding at least one more ciliurn to a 2-point tree. Considering the unique path in 
the tree between the two cilia, we see that the new ciliurn opens at least two new faces, that minimum being 
realized when the path is made of bicolored edges all of the same type. Therefore every higher-than-2-point tree 
map scales at most as N~ 2 . 


□ 


D. Characterizing LO maps 

In this subsection we prove the results for the vacuum case only, the 2-point case being similar. 

Definition 1. A map saturating the bound of the theorem is said to be LO (which stands for leading order ). For 
simplicity we shall again restrict to vacuum or 2-point maps and define a LO connected (vacuum or 2-point) map as 
a map saturating the bound, hence with f l(M) = —4 (vacuum case) or H(M) = 0 (2-point case). Trees are LO. 
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Lemma 2. If a connected map M. is LO, all its connected submaps AT C M. must also be LO. 

Proof Suppose AT C A4 is vacuum and not LO. We can pick a spanning tree T CM whose restriction T' to AT 
is a spanning tree of AT. Then fl(A4) = f l(T) = fi(T') < H(AT). 

• Let A4 = AT U (T\ T') be the map obtained from A4 by deleting the edges which are neither in T nor in Ad'. 
Performing those deletions one after another one, one gets a sequence of maps which are all connected and by 
Lemma[lJ Sl(Ad) > fl(Ai). (Notice that here it is not necessary to pay attention to whether the deleted edges 
are monocolored or bicolored.) 

• Deleting the edges of T\T' from Ad then does not change the value of Ll. Indeed, Ad \ Ad' = T\T' is a forest of 
trees, each of which can be rooted at the vertex which meets with T' (there is only one such vertex; otherwise, 
there would be a path from one to another going through T' and another path along T\T' which is impossible 
as T' C T). It is then easy to check that deleting a leaf together with its incident edge leaves f 1 invariant. By 
induction from the leaves to the roots of T \ T', it comes that fl(Ad) > fl(Ad) = fl(Ad'). 

□ 

We now want to characterize the LO maps. The above proposition implies that: 

Proposition 1. For any connected LO map, 

1. Any connected submap made of monocolored edges is a tree. 

2. For all i ^ 1, any connected submap formed by bicolored edges with colors 1 i is planar. 

3. The connected submaps made of bicolored edges only are planar. 

Proof By Lemma [2] it is sufficient to prove that the connected LO maps with only monocolored edges are trees, 
and that the LO maps with only bicolored edges are planar. 

• Item [I] The LO maps formed by monocolored edges correspond to the melonic sector of the quartic tensor 
models and consist of trees in the IF-representation RESj . 

• Item [I| It is a direct consequence of Lemma [2] applied to Item [3] of the proposition. However, we can offer a 
simple and direct proof, so that we can use this result to prove Item [3] We set * = 2 for definiteness, and denote 
A! 12 the connected map. As an ordinary map (that is, forgetting the edge colors), we denote the number of 
faces, edges and vertices F 12 , E 42 and V 12 . Notice that its faces with colors 01 and those with colors 02 coincide, 
which implies F 01 = T 02 = -F 12 • Moreover, each vertex contributes to one face with colors 03 and one with colors 
04. Therefore F 03 = F 04 = V 12 . Consequently, the function fl evaluated on such a connected submap reduces to 

H(A4l2) = 2 E 12 — Fq 1 — Fq2 ~ F 03 — Fq 4 = —2(Fi2 — E\2 + V 12 ) = 4^12 — 4, 

where the genus g\ 2 > 0 of A4i2 has been introduced. Therefore, minimizing H(A4i2) is equivalent to the map 
being planar, g\ 2 = 0. 

• Item [3] Denote A4 a connected map made of bicolored edges only. When seen as an ordinary combinatorial 
map, M has V vertices, E = Y^j =2 j edges and F faces and we want to show that its genus g vanishes 
whenever A4 is a LO map. Notice that the F faces coincide with the F 01 faces of colors 01. Therefore Euler’s 
formula for the genus can be used to write the evaluation of II on A4, 

4 

H(A4) = 2E - F - F 02 - F 03 - F 04 = 2g - 2 + V + Y,( E 0 - F o)‘ ( 17 ) 

J=2 

Then, let j £ {2,3,4} and consider the submap Mj made of edges of color type 1 j. From the above Item [ 2 ] 
for each connected component of A4 j Euler’s formula holds with trivial genus. We denote pj the number of 
connected components. A vertex of A4 which is not incident to an edge of type lj is counted as a connected 
component with vanishing genus - it has a single vertex and a single face (we have to count those isolated 
vertices since such a vertex carries a face of color 0 j and these are the ones we want to count). Therefore, the 
total number of vertices of A4 j is V and 


4 

=> Q(M) = 2g- 2 + 2(2V- E pj)- 

3= 2 


Eij — Fqj — —2 pj + V 


( 18 ) 
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Finally we show that YPj =2 Pj — 2T +1- This is proved by induction on V. It is true for V = 1 since then pj = 1 
for j = 2,3,4. Assume the bound holds for any connected maps over V vertices and let Ad be a connected map 
over V + 1 vertices. Let v be an arbitrary vertex of Ad and let AT be the map obtained by removing v and 
all its incident edges from Ad. We denote pm = Y j =2 Pj f° r Td and PM' the same quantity for AT. There is 
at least one edge connecting v to AT. If there is exactly one, say of type 12, then v (and its possible incident 
loop edges) counts as a connected component for AT 3 and AT 4 , and pm = Pm 1 + 2. If there are more than 
one, then it might be that v belongs to connected components of several color types which intersect AT, which 
implies pm < PM' + 2. The induction hypothesis on Ad', i.e. pm 1 < 2T + 1, then leads to pm < 2(V + 1) + 1, 
as desired. 

Therefore, the bound ]Cj =2 Pj — 2T + 1 implies that 

(19) 

from which is it concluded that a map of genus g > 1 cannot be a LO map. 

□ 

We also want to analyze how the various monocolored and bicolored connected submaps are connected to one 
another in a LO map. First, we show a stronger statement than just saying that monocolored edges form forests. 

Proposition 2. All monocolored edges of a LO map are cut-edges. 

This can be seen as a refinement of Lemma [ 2 ] in the case where the submap Ad' is formed of monocolored edges, 
for which the strict bound of the case [T] of Lemma |T] applies. 

Proof Let Ad be a (connected) LO map and assume that e is monocolored and not a cut-edge. Then its deletion 
results in a map Ad' which remains connected. Moreover, from the case [l] of Lemma [lj Sd(Ad') < If (Ad). This is a 
contradiction. □ 


1. LO Maps, Restricted quartic model 


In the restricted case, the maps are made of monocolored edges and bicolored edges of color type 12 only. We 
introduce some notations and definitions. 

Definition 2. Let Ad be a vacuum, connected map. 

• Removing all monocolored edges, we get a map whose connected components are submaps made of bicolored 
edges only (including the trivial submap made of an isolated vertex). Let Adi 2 be the set of those connected 
components. 

• Denote £ the set of monocolored edges. Notice that the vertices incident to a monocolored edge belong to at least 
one element of Ni 12 , by definition. 

• Let Gm be the graph whose vertices are the elements of Ad 12 and whose edge set is £. 

Theorem 2. The vacuum, connected LO maps are the maps Ad such that Ad 12 consists of planar maps of color type 
12 and whose monocolored edges are cut-edges. Moreover, the graph Gm is a tree. 

Proof We already know that if Ad is LO, then Ad 12 consists of planar components and monocolored edges are 
cut-edges. Therefore, we only have to prove that those two constraints are sufficient. Let Ad be such a map. 

Let e be a monocolored edge in Ad. It gives rise to a unique edge in Gm and the other way around. Therefore we 
identify both edges, denoted e. The fact that it is a cut-edge in Ad implies that it also is a cut-edge in Gm- Since 
this holds for all edges of Gm > it means that it is a tree. 

It remains to show that the degree of Ad is fl(Ad) = —4. We will use the tree Gm to that purpose. Denote 
Ad 12 = {Adj^lpe-R where I? is a finite set, and every Ad^ is a connected bicolored (planar) submap of Ad as well 
as a vertex of Gm- We equip Gm with a root vertex, i.e. a distinguished p* £ R and a corresponding root submap 
That induces a partial order relation on Ad 12 . m[ p P is declared the largest element and its adjacent vertices 
are all smaller, and so on down to the leaves. 

The leaves are the elements which have no elements smaller than them. In Ad, they correspond to submaps incident 
to a single monocolored edge. For every one of them, this edge is separating. Consequently, the last part of Lemma 
[T] applies, and 


D(Ad (1) ) = fl(Ad) - fi(Ad^) - 4. 


(20) 
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Here Ad^ is the submap corresponding to the leaf of Gm under consideration, which is incident to the edge e. Ad^ 1 -* 
is the map obtained by deleting Ad^ and e from Ad. Since Ad^ is planar, = — 4, hence 

= Q(M). 

The graph G M (i) is precisely G_m with the leaf and its incident edge removed. Thus, it also is a rooted tree. One 
then proceeds inductively, from the leaves to the root of Gm ■ This inductively removes all the connected pieces Ad \ f '}, 
p £ R, and their planarity is crucial so that the degree id does not increase along the removal process. 

The induction reduces Gm to a single vertex (the root) which corresponds in Ad to either a single vertex or a 
non-trivial connected planar submap of color type 12. The degree is —4 and this is the degree of Ad too. □ 

Remark. An alternative to the statement “Gm is a tree” (and its proof) of Theorem [2] consists in the following. 
Let 7~m C Ad be the submap whose vertices are the vertices of Ad, which contains all monocolored edges of Ad as well 
as the bicolored edges of a spanning tree of every maximally connected submaps Ad^ of color type 12 (not reduced 
to a single vertex). Then Tm is a tree if and only if Ad is a LO map. 

Indeed, first assume that Ad is LO and that Tm has a cycle. If it contains at least one monocolored edge, Ad cannot 
be LO. Therefore it must be a cycle whose edges are bicolored of color type 12. Obviously they belong to a single 
connected component Ad^ for some p £ R. This is however impossible since Tm only contains a spanning tree of 
Ad y). As a consequence, Tm is a tree if Ad is LO. 

Furthermore, assume that Tm is a tree (this implies that monocolored edges are cut-edges) and that Ad^ is planar 
for all p. By construction, the difference between Tm and Ad itself is a set of bicolored edges. Every one of them 
belongs to a unique maximally connected component of color type 12. Let e G Ad ^ be such an edge and let T \2 C Tm 
be the spanning tree of Ad ^ • Since the full component is planar and since e is not a cut-edge, removing it does not 
change the value of on Ad^\ which eventually gives ^(Ad^) = O {Tyf *) for all p. Moreover, the components of 
{Ad^} are connected together by monocolored edges which are cut-edges. Therefore the last statement of Lemma [l] 
applies and easily leads to fl(Ad) = LI(Tm)- As we know from the proof of the l/N expansion, the degree H is —4 for 
trees, hence fl(Ad) = —4. A similar argument will be used in the full quartic model below. 


2. LO Maps, Full quartic model 


Consider a vacuum LO map Ad in the full quartic model. We know that monocolored edges are cut-edges and that 
the bicolored components are planar. It remains to characterize the way the bicolored planar components can be 
attached to one another. Indeed, in contrast with the restricted quartic case, there are three types of bicolored edges 
and not all planar maps with bicolored edges are LO. 


Definition 3. Let Ad be a vacuum, connected map. 

• Remove all edges but those of color type 12 (respectively 13 and 14) as well as the isolated vertices this cre¬ 
ates, and denote the connected components thus obtained Ad^,..., (respectively Ad ^,..., and 

Ad«...,Adf 4 4) j. 

• Pick up a spanning tree T ^ for each i = 2,3,4 and pi = 1,..., Let Tm C Ad be the submap which contains 
all those spanning trees as well as all monocolored edges. 

Proposition 3. Let Ad be a vacuum, connected, LO map of the full quartic model. Then Tm a tree. 


Proof Tm obviously contains all the vertices of Ad and is a connected map. Let us assume that Ad is LO while 
Tm has a cycle. This cycle cannot contain a monocolored edge. It is therefore assumed to be made of bicolored edges 
only. Furthermore, those bicolored edges cannot be all of the same color type. Indeed, if that were the case, then 
they would all belong to a single connected component Ad ( j(-'' 1 (by definition of the latter and because the cycle is 
connected), but this is impossible as only a spanning tree of M[ p ^ is part of Tm- 

The cycle, say i , thus has at least two bicolored edges of two different color types. As a vacuum map itself, it is 
easy to check that its degree takes a non-LO value, Ll(£) > —2. According to Lemma |2j it cannot be a submap of a 
LO map, which is a contradiction. □ 

Proposition [3] puts explicit restrictions on the gluing of the planar components in the large N limit. In particular, 
Ad^ 1 -* and Ady ^ can at most share one vertex. If not, there would be a path in T}['^ and another path in T^ 3 ^ 
joining the same two vertices, thus creating a cycle in Tm- 
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FIG. 4. This represents the structure of the LO maps in the full quartic model. The grey areas are connected components of 
given color types. A bicolored connected component can be attached to another one on a single vertex, without forming cycles 
of such components. 


Furthermore, there cannot be “closed chains” of planar components of different color types. If and are 

not incident to one another, then they can be connected through either monocolored (cut-)edges, or other connected 
submaps M whose removal disconnects M. 

Figure [4] therefore shows the structure of the LO maps. They are planar, and made of trees of monocolored edges 
which connect bicolored connected objects. The latter can touch one another at a single vertex at most and do not 
form closed chains, thus displaying a “cactus” structure. 


IV. TREES OF NECKLACES AND TREES OF DISKS 
A. Trees of necklaces as generalizations of the quartic interactions 

Results established in the frame of a model with fixed interactions are expected to hold beyond that particular 
model. This is the content of the expected universality (in the sense used in statistical mechanics): changing the 
microscopic details (the form of the building blocks) does not affect the critical properties. This is quite well verified 
in matrix models. There, it is intuitive that moving from triangulations to quadrangulations or hexangulations is 
imperceptible in the scaling limit (mixing various interactions also brings additional degrees of freedom which upon 
fine-tuning leads to multicritical behaviors). We know this is also true in tensor models equipped with their standard 
scaling. 

Even off-criticality, one can draw interesting conclusions about generic models from the study of a single one. In 
the 2D case, for instance, this is based on the observation that hexagons can be obtained by merging two squares. 
We will apply this technique to the restricted case so as to extend the range of applicability of our results to a family 
of models. The restricted case indeed has a quite natural generalization: 

• change the quartic melonic interactions into arbitrary melonic interactions 

Vmeionic = exp-N 3 ^ X B B(T, T), (21) 

melonic B 


• change the quartic necklace interaction into a finite linear superposition of necklaces of the same color type but 
of arbitrary length 


-£>p n 2+ pb[i\ T,T). 


^necklace( 12) GXp 


( 22 ) 
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The bubble corresponding to B[^ has 2 p vertices and a single cycle made of edges alternating the colors 1 and 
3. Moreover, for each edge of color 1, there is an edge of color 2 between the same vertices. Similarly, for each 
edge of color 3 of the cycle, there is an edge of color 4 between the same two vertices. Note that for p = 1 and 
2, this reproduces the Gaussian and necklace term of the measure dp res . 

If one turns off the necklace couplings (p p ), this is just the ordinary tensor models with melonic interactions (which 
we know how to solve). If the melonic couplings are turned off on the other hand, this reduces to an ordinary matrix 
model. Indeed, defining Mab = VN T ni „ 2 „ 3n4 whose matrix indices are pairs of tensor indices, A = ( 711 , 712 ) and 
B = ( 713 , 714 ), it is found that the necklace polynomial of degree 71 in T is 

N 2+p b[ p 2 \ T.T) = N 2 Tr(MM t ) p . 

Finally, we can embed the above families of melons and necklaces into a larger one, for which we give a recursive 
definition. 

Definition 4. We say that a necklace (of color type 12) is open on the color i, for i € {1, 2,3,4}, if an edge of color i 
has been cut to form two half-edges. A tree of necklaces of type {pi,... ,p n ,p n + 1 } is obtained from a tree of necklaces 
of type {pi,... ,p n } by removing any edge of color i and replacing it with the necklace of size p n +1 open on an edge 
of color i (and preserving bipartiteness). We call this process the insertion of a necklace (see F*pwe[5|j. 

The tree of necklaces of type {p± } is simply the necklace of size p\ . The insertion of a necklace of size p = 1 simply 
is a melonic insertion. Therefore we indeed reproduce all the melonic polynomials and all the necklaces of color type 
12. Notice that the data {pi,... ,p n } does not capture the full structure of the observable. It only records the sizes 
of the necklaces which are inserted one after the other one. This will be sufficient for the enumeration of the LO 
contributions. 

Let us denote a generic tree of necklaces by C. If it is of type {pi, ... ,p n }, we set 

n n 

w(C) = ^(2 +pk) 3 (n — 1) = 3 — n + '^Pk- (23) 

k =1 k =1 

This is the enhancement trees of necklaces require to contribute at large aQ The class of model we will analyze is 
characterized by the measure 

dp(T,T)=exp(-^7V-( £ )f £ I? £ (T,T)) dp 0 ( T,T). (24) 

c 

We recall that dp o is the Gaussian measure <©• The sum in the exponential is over a finite number of trees of 
necklaces. 


Theorem 3. Let us denote the expectation of the necklace of size p as 


C p = 


N 2 +p 

N 4 




N 2+p J dp(T,T)B[l\T.T) 
^V 3 Jdp( T,T) 


(25) 


7 Notice that for each k, 2+p*. is the exponent of N expected for the necklace of size p The reason for the term —3 (n— 1) will appear in 
the proof of Lemma[3] Let just say that those observables will be represented as trees within the quartic model with (n— 1) monocolored 
edges, each of them being supposed to come with a factor N~ 3 . 
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Then the expectation of any tree of necklaces 

]y“(^(Pl.Pn}) 

]y4 \Mpi 


factorizes in the large N limit like 

n 

,...^}(T,T)> = n C Pk . 

k =1 


(26) 


The strategy to prove this theorem will be to reduce it to the quartic case defined by the measure (12), for which 
the proof is simple (given the results of the section III). To do so, we need to introduce the notion of boundary graph. 


B. Reduction to the quartic case 


The reduction from a generic model to the quartic model (with the standard scaling of tensor models) has already 
been used in Here we simply extend it to enhanced trees of necklaces. 

An open Feynman graph is defined here as a connected, bipartite, edge-colored graph where all vertices have incident 
edges of colors 1, 2, 3, 4, together with some edges of color 0 which connect some black to some white vertices. In 
particular, some vertices are not incident to the color 0. We call them free vertices. 

A Feynman graph which is open has closed, internal faces as well as broken faces with alternating colors 0, i, which 
go from a white free vertex to a black free vertex. 

Definition 5. The boundary graph dQ of an open Feynman graph Q with 2 p free vertices is defined as follows. Each 
(black or white) free vertex h gives rise to a (black or white) vertex Vh■ An edge of color i is drawn between two 
vertices Vh,v g in dQ if there is a broken face in Q which starts and ends at the free vertices h and g. The boundary 
graph is bipartite, regular of degree f, with colored edge and has 2 p vertices, but not necessarily connected. 

Note that a graph whose vertices are all free is its own boundary graph. The connected boundary graphs are 
precisely the observables of tensor models, i.e. the bubbles representing the connected tensor invariants. We recall 
that the latter appear as the 0-bubbles of Feynman graphs. Therefore, it makes sense graphically to trade an open, 
connected subgraph of a closed, connected Feynman graph for its boundary graph if the latter is connected: the result 
is still a connected, 5-regular, edge-colored bipartite graph. Furthermore, by definition of the boundary graph, this 
operation preserves all the faces which go through the subgraph but are not restricted to it (i.e. the broken faces of 
the subgraph). 


Lemma 3. Let T and T r 


be the sets of vacuum Feynman graphs generated by respectively the model defined by dp 
(12|. There is a map Q which to every graph Q £ T associates a graph Q(Q) £ F res 


(24) and that defined by dg 
such that the degrees 01(G) and Ol(Q(Q)) (calculated in the appropriate models) coincide. 


Proof We first build the map Q in a purely graphical way. This is done by exhibiting a family of open Feynman 
graphs of the quartic model whose boundary graphs are the trees of necklaces. 

Let us start with a single necklace of size p. It can be obtained as the boundary graph of a loop of p quartic 
necklaces of color type 12 connected by edges of color 0 (drawn as dashed edges), 


3 


1 





(27) 


The loop creates one closed face of colors 01 and another of colors 02 and leaves 2 p free vertices. In the boundary 
graph, each quartic necklace gives rise to an edge of color 1 and an edge of color 2. The edges of colors 3 and 4 
correspond to broken faces shared by couples of quartic necklaces. 

Assume we have constructed an open Feynman graph Qc for every tree of necklaces with n insertions, such that 
dQc = £■ We consider a tree of necklaces L n +\ with one more insertion. It comes from inserting an open necklace of 
size Pn+ i on an edge e of color i of a tree of necklaces C n of type {pi,... ,p n }- Consider its associated open Feynman 
graph Qc n • The edge e is represented in Q by a broken face of colors 0 i. The white vertex to which e is incident is 
represented in Qc n by a white free vertex (no edge of color 0 incident). Let us attach to this free vertex the following 
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piece of graph 



(28) 


where the loop of quartic necklaces contains p n +\ of them. Denote Gc n +i the open graph obtained this way. We 
want to show that dGc„+i = Ai+i- Since the edges ei, e 2 of color 0 are cut-edges, one can replace the subgraph Qc n 
with its boundary graph C n and also the loop of quartic necklaces in (281 with its boundary graph (a necklace of size 
Pn+i)- Denote Q n +\ the graph obtained this way, then dGc n+1 = dQ n + 1 - The only portion of Gn+i which contains 
non-free vertices looks like 



(29) 


By applying the definition of the boundary graph (and separating the cases i = 1,2 and i = 3,4), it is easy to see 
that it creates a necklace insertion of size p n +i onto e in C n . as desired. 

As a consequence, given a graph 5eJ, one can replace all its trees of necklaces {£} with their associated open 
graphs {Gc}- This way, one obtains a vacuum Feynman graph Q(G) € F res from the restricted quartic model. 

If G £ -A contains £(G) edges of color 0, F oi (Q), i S {1,2,3,4}, faces of colors Oi, and a set of trees of necklaces {£}, 
it comes with weight where 


- fi(0) = E F 0l (Q) - 3 £(Q) + E u(£). 

c<zg 


(30) 


i=1 


From the definition of uj(jC) in (23), 

4 n(C) 

-n(G) = E F 0 i($) - + EE < 2 + p *09) - 3("(9 - x ) 

i =1 CCQ k—1 

n(C) 4 p n(£) 

= -s [m + e ( 2 (n(c) -1) + E Pk (£))) +E F oi(G) + E 3 H9 - 1) + E ( 2 + 4 P*0C)) 

£cS fc=1 j=l £cs L fc=1 

Now we partition the above contributions in order to assign them naturally to Q(G). We recognize 

n(C) 


(31) 


£(Q(G)) = £(G) + E ( 2 «(9 - 2 + E Pk{£)) 


(32) 


c eg 


fc=i 


as the number of edges of color 0 of Q(G): it still has those of G, plus those of the graphs Gc- In one of those graphs, 
we count pk edges of color 0 to form each loop of quartic necklaces whose boundary graph is a necklace of size pk ■ If it 
has n(£) necklaces, there are 2 n(£) — 2 edges of color 0 added to connect them (like e\, e 2 in (28)). We also recognize 


n(£) 

bi 2 (Q(G)) = E E p*@)> and kq@)) = E «(9 - !> 

c<zg k =i c<zg 

the number of quartic necklaces and the number of melonic quartic bubbles of Q(G)- Finally 

Foi(Q(G)) = Foi(G) + (Si ,i + Si, 2 ) E n (9i 

ccg 


(33) 


(34) 
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is the number of faces of Q(G) of colors 0 i, since every necklace represented as in (28) has two internal faces, one of 
colors 01, and one of colors 02. 

It comes 


4 

- n(g) = -M(Q(g)) + 4 MQ(S)) + 3&(Q(S)) + F ^Q^)) = -fi(Q(S)), (35) 

i =1 

as assigned in the restricted quartic model to Q(G). □ 


C. Factorization 


Proof of Theorem [3] The expectation of an observable is a sum over all the closed Feynman graphs generated by 
the model which contain the observable as a marked subgraph. We choose the observable to be a tree of necklaces. The 
graphs contributing to its expectation are closed graphs made of trees of necklaces connected through edges of color 
0. We apply the map Q to represent them in the quartic model and then use the intermediate-field representation. 
Through this process, a tree of necklace £ of type {pi,... ,p n } becomes an open graph Gc which in the IF representation 
becomes a map Tc with n bicolored loops of sizes pi ,..., p n connected in a tree-like fashion (no closed circuits between 
those loops) via n — 1 monocolored edges. 

From the combinatorial description of the LO graphs of the restricted quartic model, it is clear that there exist 
graphs Q contributing to the expectation such that Q(G) is a LO graph. From Lemma [3j it thus comes that Q 
contributes to the LO of the expectation if and only if Q(G) is LO in the quartic model. 

Let G be a LO graph contributing to the expectation of £ which is of type {pi,... ,p n }. In the IF representation, 
Q(G) contains Gc as a marked subgraph with n — 1 monocolored edges and n bicolored loops. All those monocolored 
edges are cut-edges in Q{G). Cutting them yields n connected components, each one containing one of the n bicolored 
loops as a marked subgraph. This establishes a factorization of the expectation of £ onto its necklaces, 

_/V w (£{pi._ AC _ 

- w -<A Pl ,...,Pn}(T,T)> = l[C Pk . (36) 

k —1 


Moreover, any set of graphs contributing to the expectations of the loops of sizes Pk {£), k = 1,... ,n, gives rise to a 
LO contribution to the expectation of £ by gluing back the monocolored edges. This proves that the contribution of 
a necklace is C Pk = C Pk as defined in Theorem [3j □ 

Although this result is simpler to prove using the IF representation, it can also be understood directly in terms of 
the edge-colored graphs generated by the initial model. If £ is the observable whose expectation we are interested 
in, let £ Pk be one of its necklaces, and let v,v € £ Pk be white and black vertices incident to edges of color i which 
connect it to other necklaces. By construction of £, v and v can be chosen such that their incident edges of color i 
are 2-cut-edges. In G, there are also two edges of color 0 incident to v and v. 


The factorization (36) means that any path which starts at v (or v) and follows its incident edge of color 0 does 
not hit £ before it comes back to a vertex of £ r . 


-Pk ■ 


• Moreover, the equality C Pk = C Pk can be found by studying the number of faces of colors 0* which go through 
£ Pk . The unique face of colors 0* which crosses v follows the edge of color i incident to v and necessarily goes 
along the edge of color i incident to v. It happens this way because the two edges of color i incident to v and v 
are 2-cut-edges, i.e. they carry a 2-point function with external edges of color i. Therefore, when counting the 
number of faces 0 i which go through £ Pk , everything happens as if v and v were connected by an edge of color 
i. This means that the counting can be performed directly on the isolated necklace of size pu, which yields the 
desired equality. 


D. Enumeration: Large N free energy and expectations 

Thanks to the combinatorial analysis of the quartic model performed in Section m the problem of evaluating the 
expectations of trees of necklaces at large N has been reduced to determining the expectations of the necklaces: this 
is the content of Theorem [3j Calculating those expectations amount to enumerating the number of Feynman graphs 
containing a necklace of fixed size as a marked subgraph. To do so, we will write an equation a la Tutte, known in 
our context as Schwinger-Dyson equations. 
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The Schwinger-Dyson equations of random tensor models have been described in 1491 . and the structure and large 
N solution have been described in 150 El]. We recall that an observable is a U (JV) 4 -invariant polynomial which can 
be represented as a connected, bipartite, 4-colored graph. Then define the following two operations. 

• The contraction of an observable B with 2 p vertices along the black and white vertices v,v £ B, denoted B / ( v , v) 
is obtained by connecting v and v with an edge of color 0 and taking the boundary graph. Equivalently, one 
removes v and v and reconnects the edges while respecting the colors, to get a set of disjoint observables on 
2 (p — 1 ) vertices. 

• The composition of two observables B,B' along the pair v £ B,v' £ B' is obtained by connecting v to v' via 
an edge of color 0 and then taking the boundary graph. Equivalently, one removes v from B and v' from B' 
and reconnects the edges while respecting the colors. If B has 2p vertices and B' has 2 p' vertices, then their 
composition, denoted B *(«/,„) B', has 2(p + p' — 1) vertices. 

Let B be an observable, represented as a connected, bipartite, 4-colored graph, with a marked white vertex v. There 
is one Schwinger-Dyson equation for each such observable with a marked vertex. It reads 


Y v)) - N 3 (B) - J2 N U(C) t c Y, *(®M0 B) = 0, 

veB c v'ec 


where {£} denotes the set of trees of necklaces which defines the measure (24). 

We choose B = B[^ the necklace of size p, and multiply the equation by Ni p ~ 4 ) /N 4 . Then 


, , ]\f2+(p-k-l ) ^ 2 +k 

-{B$/(v,v)) = - w ~^(B[ 2 } ) = CV fc _rC fc , 


• If v and v are separated by k white vertices, then 

Nip- 1 ) 

N 4 

which is the same large N factorization as in ordinary matrix models. 

• The middle term of the equation, which comes from the Gaussian part of the measure, reads 

N^- 1 ) , x 

iV 3 ^) = C p . 


N 4 


(37) 


(38) 


(39) 


• If £ is an interaction of the model, of type {pi,..., p n } and v' £ C belongs in a necklace C p ., then the composition 
with B [ P 2 increases its size to Pj + p — 1 , and 

Ni p — 1)+ W (^'{P1.Pn>) _/y“(£{p 1.Pj+P-l,...,Pn}) 

(B*(v',v) £{pi,...,p„}) = —.Pi+P-I, —,Pn>) 


= c, 


n 


pj+p-i ii w p/=> 

k=l,...,n 


(40) 


where the last equality makes use of Theorem [3j 


Gathering the above pieces, equations (371 become a system which determines the rescaled expectations (Ci, C 2 , C 3 ,...). 


It takes the following generic form. Let^ v\xi, x%, X 3 ,...) be a polynomial in a finite number of variables (aq, X 2 , £ 3 ,...), 
and denote d p the derivative with respect to x p . Then the Schwinger-Dyson equation has the form 


p - 1 

C P = Y Ck C P ~ k ~ 1 + Y. 3 djV{C u C 2 , C 3 ,...) Cj+p-i. 


(41) 


k =0 


j> 1 


This equation has a clear interpretation in terms of branching planar maps (i.e. trees of disks). C p is the number 
of maps in the class with root vertex of degree p, or by duality, boundary face of degree p. The quadratic term 
corresponds, as usual for equations a la Tutte , to the case where the root edge is a bridge. If not, its removal extends 
the length of the boundary face from p to p + j — 1, which is also usual for planar maps. However, it here comes with a 
weight jdjV(Ci,C 2 , ■ ■ ■)■ When this is independent of the numbers { Ck} (the linear terms of V), one simply recovers 
Tutte’s equation for planar maps with bounded face degrees. All other terms of V correspond to gluing “on top” 
of our map other maps with prescribed boundary face degrees. This superimposition generates a branching process 


which we already analyzed at the combinatorial level in the quartic case in Section III in terms of trees of disks. 












18 


Equation (41) can be found in the existing literature devoted to multi-trace matrix models 


Therefore, 


we will only describe its critical properties without reproducing the detailed analysis. The free energy behaves like 
(g — g c ) 2 ^, where g is the “cosmological constant” (an overall coupling, or a coupling of ordinary single-trace terms 
of the potential), g c is the radius of convergence of the generating function for (C p ) p >i (it depends on all the other 
couplings of the model) and 7 is the entropy exponent. The latter characterizes the proliferation of some combinatorial 
degrees of freedom of the maps like n 7-3 where n is the size of the objects. It classifies the various universality classes 
which can be achieved. 


• When the planar components become critical while the branching process stays off-criticality, the latter is washed 
away in the scaling limit and the string susceptibility exponent of pure 2D gravity is recovered, 7 = —1/2. 

• If the situation is reversed (non-critical maps, critical branching), 7 = 1/2, as expected for trees since the planar 
components become irrelevant in the scaling limit; this is the universality class of branched polymers. 

• If both are simultaneously critical (a phase reached by tuning one additional coupling to a critical value), the 
exponent 7 = 1/3 can be reached. This phase transition describes a proliferation of baby universes. 

• Tuning more couplings leads to more phases, with exponents 7 = p/(n + m+ 1), where p < n and m are integers. 
It is argued in 03] that for n > 1 the polymerization is too strong for the surface to support macroscopic loops. 


V. ANALYTICITY OF THE MODEL 

We discuss the analyticity of the 2-point function G 2 of the general quartic model in a cardioid domain of the 
complex plane, using the Loop Vertex Expansion. The proof, being similar to the one of the non-enhanced case m , 
will be discussed in less details. 

Theorem 4. The 2-point function Gi of the maximally rescaled general quartic model with coupling constants A c = 
\Xc\e l ^ c is analytic for |Ac| < yg cos 2 (f> = maxe(</>c). 

To prove this theorem, we will first use a new series expansion of log Z in the intermediate field representation, (the 
Loop Vertex Expansion), then prove the absolute convergence of the series for all A’s inside the cardioid domain. 


A. Loop vertex expansion 


Using a standard Hubbard-Stratonovich transformation on (11), the model can be re-written in terms of seven 
hermitian matrices a c , where for i £ {1,..., 4}, a 1 is a N x N matrix, and for i > 2, a 11 is a four indices tensor seen 
as an hermitian N 2 x N 2 matrix acting on pairs of indices of color {1, i}. The measure on those matrices (known as 
intermediate fields) is 


dp,(cr) = dgo(cr)e 


-Tr^InJV’ + Etl V^( l2,U °® CTl ) + 


(42) 


where l fa is the identity matrix over the set of indices B and dpo(a) is the normalized Gaussian measure of covariance 
1 over the a matrices. 

Denoting A(a) = *Ec \J (g) cr c ) and R(a) = [l D + A\ \ the generating function of the cumulants 
(fl3| is 


Z(J, J) = / e 


- p N- 3 T. n .n^R(cy)nnJr 


dp,(a). 


(43) 


Hence the 2-point function can be expressed as the connected expectation value of the resolvent, i.e. 

G » = | < 44 > 
Moreover, it admits an expansion as a sum of trees, via the BKAR formula [581 159| . For any tree T v with v vertices, 

• to any edge £ of T v , we assign a parameter ut £ [0,1], 

• for any pair of vertices {*,/}, P%j being the unique path joining the vertices, we define w i: j = min ,eePi u e, and 
we define wu = 1 , 
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• we define a new set of a matrices {a 1, c } for each vertex i of T v , 

• we define the interpolated Gaussian measure Ht v ,u as f cr^n a mm, dpT v ,u(&) = Wijdnmdnm- 

Then we can express G 2 as a sum over plane trees T v tca)},i* with v vertices, colored edges t and a single ciliated 
vertex i*. Each corner p of a vertex i bears a resolvent R{cr l ) whose indices , nf p (c £ {1..4}) are contracted along 

the faces of color c of the tree. Therefore, the vertices have a structure of 4-stranded vertices and to an edge £ij 
between i and j corresponds the contraction 


_ Aj c c \ j § (S -D\C T>\C ] . 

\ n i,q Tl i,q + lJ '<1 D,P+1 J,P *,«+1 \ rl j, p Tl j,p+lJ 


(45) 


The ciliated vertex bears a cilium at a corner, that, for the 2-point function, is just a trivial contraction of the 
adjacent strands For higher order cumulants, the strands have to be re-connected according to non¬ 

trivial permutations over cilia, and with a factor involving Weingarten functions HQ HU® [EE]. The cilium acts as 
an edge for the ciliated vertex i*. increasing its degree and number of corners by one. Thus, 


<&)«. = ts£ 1 


£ 


N 4 ^ vl 

v > 1 Oj,{C( 0 },i* 


][J dun 


leT v 


degree(i) 

dfj-TpA^Yi ( n n 

n,n \i£V(T„) P= 1 


x ( TT 2Ac W 
X II AT4-ICWI d 

/ 


(46) 



FIG. 6. A 4-stranded tree with a cilium. Resolvents are represented as thick bars crossing each strand of a corner. 


B. Cardioid domain 


The series (46) is absolutely convergent for coupling constants |A c |, \(j)c\ < 7r inside a cardioid domain |Ac| < 
p cos 2 ( 2 )’ ^ = max C 0C- This can be shown by bounding the contribution of each tree 70 i {c(^)},i* in (46) using 
iterated Cauchy-Schwarz inequalities m ■ Writing 




d(») 

n n«(A 

uev(T„)p =1 


n 


(47) 




a colored plane tree T made of 4-stranded vertices but with no resolvent is structurally an intermediate field Feynman 
graph and thus comes with the amplitude A(T) as in (15). Resolvents can be taken care of by applying a by now 


rather standard technique H3E2HS3. Starting at the cilium, we order the (2v — 1) resolvents along the clockwise 
contour walk around the tree, indexing them as Ri, R- 2 - .. Then, we draw a line along the unique path between R\ and 
R v . Everything left (respectively right) of this path is called A (resp. B), then we use the following Cauchy-Schwarz 
bound: 


C(T vdcW} ,i*) = (A\ R l ® R v ® | B) < ll-Rillll^HV {A\A) (B\B). (48) 


Within the cardioid, ||i?(cr)|| < l/cos(^/2) with $ = maxc fic- Moreover, the scalar product graphs {A\A) and 
(B\B) have the same structure of plane trees with resolvents as the original graph, except that (A\A) now bears 2v — 4 
resolvents and (B\B) has 2v — 2. Repeating the same process on {A\A) (but using arbitrary corner as i?i) will give 















20 



FIG. 7. First Cauchy-Schwarz iteration applied on the tree from Fig|6] The cilia could have been deleted as they are merely 
identity operators. 


rise to graphs bearing 2v — 6 resolvents and so on. Repeating the process until no resolvents remains on any graph, 
we have 


2v ~ 1 Ft 

c(T v , {CW} ,i*) < n iiw^in^^n^- 

i=l 0a 0 b 

/ i \ -(2v-l) 

< (cos ^ N f ^ R \ 


f (b B ) 

2 v 


(49) 


where F{Q) is the number of faces of Q, and qb are the graphs bearing no resolvents created from the right amount 
of iterations of the Cauchy-Schwarz bound on the graphs {A\A) and ( B\B) and Tv is the tree with all resolvents 
removed. The second inequality arises from the conservation of the number of faces during the Cauchy-Schwarz 
process. Indeed, at each iteration, the number of faces is multiplied by 2. 

For any tree, F(T} R ) = 4v — ^ f C(l?), therefore, 


|(G 2 )„n| <]T^cos^ n 

v>l 77,{c(«)},i* 


2 | A C (£) [ 

cos 2 | 


(50) 


Finally, using Y2t {C(t)} ■ . 1 = F 1 an< ^ ^ 2 % — Pi we find absolute convergence of the series in the cardioid 

domain with p = A 

The perturbative expansion in A can be recovered by applying Taylor expansion with integral remainder individually 
on each tree, up to an uniform order. This corresponds to expanding beyond trees by adding additional loop-edges 
to the graph. This mixed expansion [48] here does not have such simple properties as in [47], because additional 
two-colored loops do not always lead to higher orders in 1/N. The integral remainders however, still give a bound on 
the overall Taylor remainder of G 2 which is good enough to prove Borel summability of the Feynman graph expansion 
uniformly in N. 

This also extends to higher-order cumulants, introducing trees with more than one ciliurn and Weingarten functions 
[47l 118] . As in [47], recovering proper scalings for the cumulants is done by Taylor expanding further than trees before 
using iterated Cauchy-Schwarz bounds. In this case however, two-colored loop edges do not lead to automatic scaling 
improvement, thus the expansion must be done with respects to single colored couplings only. This requires introducing 
several independent coupling constants. 


VI. CONCLUSION 

Pushing further the approach of 3D] and based on the observation of |3T] that non-melonic observables can be 
enhanced, we have studied two (related) families of models with maximally-rescaled interactions. 

First we introduced the most generic quartic model at rank four, which features four melonic interactions and 
three necklaces which are matrix-like observables. In addition to the existence of the 1/iV expansion, we performed 
a combinatorial analysis of the leading order (LO) graphs in the large N limit. In the case where only a single 
necklace has a non-vanishing coupling constant (called the restricted quartic case), Theorem [ 2 ] fully characterizes the 
LO graphs. They have the structure of trees whose vertices are disks, one tree being inserted on any face of the other 
one and so on. 

The restricted case was then generalized to a generic class of interactions which we called trees of necklaces (the 
latter include all necklaces of a fixed color type as well as all melonic observables). We obtained a factorization 
theorem (Theorem [ 3 ]) which reduces the large N evaluation of the expectations of trees of necklaces. This, combined 
with Schwinger-Dyson equations, leads to the exact enumeration of the LO graphs contributing to the free energy or 
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to such expectations. A phase transition (with a family of positive entropy exponents) is found between the branched 
polymer phase and 2D quantum gravity, which can be thought of as the proliferation of baby universes. 

Beyond the specific outcome of the analysis (the Schwinger-Dyson equations reduce to the same equations as in 
matrix models with multi-trace invariants) we think the method we used is quite powerful and worth summarizing 
here. 

• It starts with a quartic model. This is because for quartic interactions the intermediate field method is straight¬ 
forward and powerfuj^] simplifying greatly the combinatorial analysis. In particular, it was first used in [55] to 
study a quartic tensor model beyond its LO and first sub-leading correction (resulting in the first double-scaling 
limit for random tensors). More recently, the intermediate field method for quartic models has been used to 
construct such models non-perturbatively. 

• Then a generic class of interactions is introduced: they are the boundary graphs of the quartic model studied 
before (in our case, of the restricted quartic model). This allows to consider the Feynman graphs of those new 
models as a sub-family of those of the quartic model. The combinatorial analysis performed in the quartic case 
thus applies. 

• That combinatorial analysis is then a fruitful input which simplifies the Schwinger-Dyson equations and can 
turn them into a solvable set of equations. 

This scheme was in fact first applied in m■ It enabled extending the double-scaling limit of the quartic melonic 
model of jj55] to arbitrary melonic models (invariant under permutations of the colors though). This is however the 
first time it is applied to non-melonic interactions. We emphasize that the Gaussian expectations of a sub-family of 
trees of necklaces were calculated recently in [55] (with a completely different method). Clearly, the present method 
goes far beyond the previous achievements. 

The present results can be continued in several directions. First, it would be interesting to know the continuous 
spaces towards which the “trees of disks” converge in the GHP sense. It is tempting to think about them as general¬ 
ization of the “looptrees” [52], but further investigation is needed. The double-scaling limit of that model could also 
be studied to understand how the double-scaling limit of ordinary tensor models can cross over to the one of ordinary 
matrix models, and how it relates to the double-scaling limit of multi-trace matrix models [451 . 

In this article, we left unsolved the full quartic model. Indeed, the presence of several matrix-like observables 
(necklaces of color types 12, 13 and 14) makes the Schwinger-Dyson equations more complicated. The fact they have 
different color types reveals that this model makes a better use of the tensorial structure. We expect it to depart from 
the multi-trace matrix models. 

Another way to further take advantage of the richness of tensorial invariants is to introduce more maximally rescaled 
models. Remark that already at rank 3 and order 6 , in addition to the six melonic interactions, there is a single non- 
melonic tensor invariant, whose bubble is the bipartite complete graph A '3 3 , which is non-planar and by Kuratowski’s 
theorem, is a kind of canonical source for non-planarity in bipartite graphs (see Figure [l]). When maximally rescaled, 
it could generate an interesting leading order. This is left to future study. 

Finally it should be interesting to define similarly rescaled models of tensor group field theories, with or without 
Boulatov type projectors, to generalize the growing list of renormalizable models studied so far [551 - f75] . In particular, 
they should allow to explore in depth the frontier between asymptotically free tensor field theories and asymptotically 
safe non-commutative field theories [HHZBI- 
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